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A method is presented for determining the stress-strain relation in shear for isotropic 
fiat sheet. The method requires measurements of deflection or of extreme fiber bending strain 
in the center portion of a twisted square plate. The octahedral stress derived from the 
stress-strain curve in shear for a twisted plate of aluminum alloy agreed within about 5 
percent with the octahedral stress derived from tensile and compressive tests of the same 
material. The difference was much greater for specimens of mild steel with a definite yield 
point. Unfortunately, the method is difficult to perform on thin sheet because of the small 
size of the specimen. An examination of alternate methods indicates particular promise for 
the use of an annulns of constant i hickness. 



1. Introduction 

There has long been an interest in the stress-strain 
relations in shear of materials in the form of sheet. 
The interest has beeome more pronounced as working 
stresses and strains have increased and as it 1ms 
become apparent that the mechanical properties of 
thin sheet could not be deduced directly from 
mechanical tests of heavier sections. The most satis- 
factory method, perhaps the only really satisfactory 
one, for determining the stress-strain curve of thin 
material in shear requires the material to be in the 
form of tubing with circular section. As long as the 
tube material is isotropic, its stress-strain curve in 
shear can be computed conveniently from the torque- 
twist curve of the tube by following the procedure 
given in fl]. 2 

There is no dearth of suggestions for methods of 
determining the stress-strain curve in shear of flat 
sheet. Figure 1 indicates 10 suggested methods for 
stressing thin sheet in shear that have come to the 
authors' attention. The first two methods deal with 
tubular specimens. 

Method (a) is nothing more than the bending of a 
strip of the sheet into a tube of large diameter, which 
is then tested like the conventional long tube for 
which a satisfactory theory exists [1]. The tube 
would have to have a large diameter so that the 
strains introduced in bending the strip would be 
small compared to the strains to be measured. 
Difficulties might be experienced with deviations 
from the theoretical uniform shearing stress caused* 
by variations in clamping forces along the edges of 
the specimen. Difficulties might be experienced also 
with the discontinuity at the joint and with the 
transmission of sufficient torque to obtain stresses 
in the plastic range. 

Method (b) was suggested to one of the authors 
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by J. M. Frankland of Chance Vought Aircraft. 
Four strips of the sheets are assembled into a tube 
of square section. The strips will be subjected to a 
uniform shearing stress by equal and opposite torques 
applied to the ends, if the torsional rigidity of the 
joints along the four corners is negligible compared 
with the torsional rigidity of the tube as a whole. 
It is probable that this corner effect is not negligible 
for tubes that are sufficiently thick-walled to trans- 
mit, without buckling, torques that extend the shear- 
ing stresses into the plastic range. 

Several methods have been suggested for stressing 
a specimen of sheet in shear by forces acting in the 
plane of the sheet. 

Method (c) is the well-known square-picture-frame 
method, in which a square specimen is held between 
four pairs of rigid clamping bars with pin joints at 
the corners. The specimen will be placed in a state 
of uniform shear by equal and opposite forces acting 
at opposite corners of the frame, provided that the 
sheet is perfectly flat, the clamping bars are infinitely 
rigid, the edge effect due to deviations from ideal 
clamping does not extend into the middle portion 
of the sheet, and the pin joints are frictionless. 

Method (d) tries to get around the difficulties with 
the edge effect in method (c) by making the speci- 
men long so that the shearing stress is practically 
uniform along the middle portion of the specimen. 
However, an end effect will remain because the shear- 
ing stress reduces to zero at the corners. This end 
effect may be evaluated by making tests of the same 
material in picture frames of different lengths. The 
procedure would be costly and time consuming. 

Method (e) was used by Coker and Filon [2] and 
Bollenrath [3] to stress thin sheet in shear. The 
loads are applied symmetrically to subject each one 
of the two "bays" of the specimen to a state of stress 
in which the shearing stress along the center line of 
each bay varies from zero at the ends to a nearly 
constant value near the middle. As in method (d), 
the end effect may be eliminated by testing speci- 
mens of different lengths. 

Method (f), suggested by W. R. Osgood, subjects 
a single strip to a state of stress which approximates 
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Figure 1. Methods suggested for stressing thin sheet in shear. 

that in each bay of the specimen for method (e). 
Method (g) avoids end effects by testing the sheet 
in the form of an "endless strip" or annulus. Equal 
and opposite torques are applied at the inner and 
outer edges of the annulus. The method has the 
advantage that, except near the periphery and the 
inner circle, the shearing stresses would be known 
accurately, for isotropic sheet, in the plastic as well 
as in the elastic range. There are some disadvan- 
tages, however. The shearing stress varies inversely 
as the square of the radius of the circle along which 
it acts. Hence there will be a nonlinear strain varia- 
tion along the gage length of any finite gage mounted 
at 45° to the tangential direction. As a result, the 
average strain along the gage will differ from the 
strain at the midpoint of the gage, unless the annulus 
is large enough to make the nonlinear variations 
negligible. If the annulus is large enough, it may be 
difficult to apply sufficient torque to stress the ma- 
terial into the plastic range. The testing machine 
would require heads that are rotatable about their 
axis only, to prevent bending of the specimen by 
rotation about other axes. The method would give 
only an average stress-strain curve for anisotropic 
sheet. 



Method (h) eliminates some of the difficulties of 
method (g). The annulus is machined into a disk 
that varies in thickness inversely as the square of 
the radius. The material in the annulus will be sub- i 
jected in that case to a uniform shearing stress by 
equal and opposite torques applied at the inner and 
outer edge. The disk must be sufficiently thick near 
the periphery to develop strains in the plastic range 
before buckling occurs and the variation in thickness 
with radius must be small enough at the radius where 
the strains are measured to assure a practically uni- 
form shearing stress throughout the thickness. A 
practical difficulty is that of machining the specimen 
accurately to the desired thickness variation with- 
out changing the properties of the material by the 
machining operation. 

The methods so far mentioned subject the sheet 
to a shearing stress that is substantially uniform 
throughout the thickness of the sheet. In the re- ' 
maining two methods the shearing stress is varied 
in the thickness direction from one extreme to an 
opposite extreme. 

In method (i) a state of shear is produced in a 
square plate by twisting the plate into an anticlastic 
surface by equal and opposite tranverse forces ap- 
plied to adjacent corners. The state of shear will 
be uniform for planes at a fixed distance from the 
median plane which will undergo zero shearing stress. 
The upper and lower faces of the plate will be under 
equal and opposite shearing stress. Objections might 
be raised against this method in so far as it neglects 
deviations from uniform shear due to various causes, 
such as having large transverse stresses near the 
corners and no shearing stresses along the edges. 
These objections will be discussed in detail in the 
body of the paper. Experimental difficulties may 
arise in applying the corner loads to the specimen 
without interfering with strain and deflection meas- 
urements on the faces of the plate. 

Method (j), suggested by A. H. Stang, uses a 
much simpler scheme for producing a state of stress 
similar to that in method (i). A specimen in the 
shape of a long strip is placed in a state of shear by 
applying equal and opposite torques to the ends of 
the strip. Shearing strains can be measured with 
strain gages or they can be obtained from the twist 
per unit length of the strip. However, if the strip is 
made long enough to eliminate end effects, it is 
likely to deflect laterally. This will superimpose 
bending stresses on the torsional stresses. The same 
objections might be raised in this case as in case (i) 
with regard to edge effects. The shearing stress at 
the extreme fiber will be practically constant in the 
center portion, far from the edges of the strip, but 
it will decrease rapidly to zero as the edges are 
approached. The value of the extreme fiber shearing 
stress near the center of the strip will depend on 
this edge effect. It can be computed exactly in the 
elastic range, but no exact solution is known for the 
plastic range. 

Consideration of the above 10 methods suggests 
that at least three of them, the square-picture-frame 
method (c), the annulus of constant thickness method 
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(g), and the twisted-square-plate method (i), should 
be Irieci out experimentally. 

An experimental study of method (c), the square 
picture frame, was conducted by Wm. R. Osgood 
about 5 years ago. Method (c) was selected because 
it appeared to have particular promise and because 
it had been used in the aircraft industry to estimate 
shear properties of sheet. A fixture had been de- 
veloped at the El Segundo Plant of the Douglas 
Aircraft Co. to make such tests. A duplicate of this 
fixture was obtained, and tests were made at the 
National Bureau of Standards for the National 
Advisory Committee for Aeronautics on several 
4.5X4.5X0.063 in. 24S-T aluminum alloy plates. 
The results were disappointing. They were reported 
to the NACA, but were not published. The strains 
were far from uniform even in the center portion of 
the plate farthest from the edges. The extreme 
fiber strains in this optimum region deviated up to 
10 percent from the average strain. The deviations 
persisted even after an attempt had been made to 
improve the conditions of clamping along the edges 
by remachining parts of the clamping fixture. 

Two years after completion of the work for the 
NACA, it was decided to try method (i), the twisted 
square plate, since Nadai had shown in 1915 [4] 
that this method provided a simple means for 
determining the shear modulus of sheet material. 
More recently the test has been extended to ortho- 
tropic materials at the Forest Products Laboratory 
[5] to determine the shear moduli of wood. 

2. Nomenclature 

The following nomenclature is used in the paper: 
a— half-span for deflection measurements. 
b = width of specimen. 
E= modulus of elasticity in tension or 

compression. 
0= modulus of elasticity in shear. 
h= thickness of specimen. 
m = moment per unit length. 
M= bending moment. 
M t = twisting moment. 
P=load at each corner of plate. 
E= radius of curvature along diagonals of 
plate. 
?/ ? r, w= displacements parallel to #-, y- and z- 
axes, respectively. 
x, ?/=axes parallel to diagonals of plate, also 
planes normal to the corresponding 
axes. 
./•', ?/ = axes parallel to edges of plate, also 
planes normal to the corresponding 
axes. 
?=axis parallel to thickness of plate, also 

plane normal to 2-axis. 
7= shearing strain. 
y = octahedral strain. 
e= normal strain. 
i> = Poisson's ratio. 
<y = normal stress. 
r— shearing stress. 
t = octahedral stress. 



3. Principle 
3.1. Elastic case 

An elastic square plate will be in a state of pure 
shear, in which the shearing stress is proportional to 
the distance from the neutral surface of the plate 
when the plate is bent into an anticlastic surface 
with equal and opposite principal curvatures. The 
bending of the plate into such an anticlastic surface 
produces equal and opposite extreme fiberbending 
stresses along the two principal axes. Nadai has 
shown that this remarkably simple state of stress 
may be realized by applying equal and opposite 
tranverse forces to adjacent corners of the plate. 

Nadai's proof [4] makes use of Saint-Venant's 
principle of the equivalence of equipollent loads. 
Nadai notes that for a load P at each corner, each 
edge of the plate is loaded by equal and opposite 
forces P/2 acting at two corners (fig. 2). The load 
along the edge will remain unchanged upon applying 
balancing forces ±P/2 at frequent intervals along 
the edge, as indicated in fig. 2. Any element of 
length of edge Ax' between two adjacent equal and 
opposite forces P/2 will then be subjected to a twist- 
ing moment (P/2) Ax/. According to Saint-Venant's 
principle, this edge load is equivalent to any edge 
load which results in a constant twisting moment 
per unit length 



1 P ^ ' P 

m ^=^2 AX= 2 J 



(1) 



uniformly distributed along the edge of the plate. 
The stress distribution, in regions of the plate suffi- 
ciently far from the edge to satisfy Saint- Venant's- 
principle, will then be that corresponding to uniform 
twisting moments -\-m x , yfj —m x , y , in the x' y y' planes, 
normal to the x' , y r axes, respectively. The twisting 
moments ±m x t y t are associated with principal bending 
moments m xi m y on the y, x planes normal to the 
diagonals of the plate. Equilibrium of moments 
acting on the edges of a triangular element of the 
plate (fig. 3) bounded by planes parallel to the x, y f 
and y' planes shows that 



m x 



-m v 



— m x 



P/2. 



(2> 



The plate is bent into an anticlastic surface by equal 
and opposite bending moments per unit length P/2 
acting on the y, x planes normal to the diagonals of 
the plate. The surface will have its planes of prin- 




Figure 2. Twisting of square plate by forces on corners* 
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Figure 3. Edge moments acting on element of twisted square 
plate, 

cipal curvature normal to the diagonals of the plate. 
If R is the radius of curvature along the diagonals, 
the deflection relative to a plane tangent to the center 
of the plate is given by 



w= 



(3) 



The distribution of strain and of stress throughout 
the plate may be computed on the assumption of the 
Kirchhoff plate theory that plane sections remain 
plane. The normal strains e z , e y must then increase 
linearly with distance from the neutral fiber 



2z 

= T €l 



2z 

= h 



€2> 



(4) 



where €i, € 2 are the principal extreme fiber strains; 
z, the distance from neutral fiber; and h, the thick- 
ness of plate. The maximum shearing strain on the 
z plane is 

2z 
V = yx'v' = ex—e y =-Y 7i> (5) 

where 71 is the shearing strain at the extreme fiber 
7i=€i— € 2 . Because of the antisymmetry of the 
deformation ci= — e 2 , so that 



7i=2 €l . 



(6) 



The shearing stress corresponding to the strain y is 
given by 

r=Gy (7) 

nside the elastic range. The twisting moment per 
unit length is given by 



Tflxlyi- 



-l 



ft/2 
ft/2 



rzdz. 



(8) 



Substituting (1), (7), (5), and (6) in this equation 
gives 



P 1 

2 = 3 €l0h * m 



(9) 



From which 



U ~2h 2 ti 



(10) 



Hence the shear modulus can be computed from the 
slope P/ei of a plot of load against extreme fiber strain 
by substituting P/ei in (10)* 

The shear modulus may be determined also from 
the deflection w a of the center of the plate relative to 
a line connecting points a distance a from the center 
along one of the diagonals, figure 2 . If plane sections 
are assumed to remain plane, the curvature 1/R along 
the diagonal x is given by the familiar formula 



1 _d 2 W_ 6i 

pr~d^~h\2: 



(11) 



where w=normal deflection. Substituting (11) in 

(9) gives 

€,a 2 3 a 2 P .. , 

W<.=— =oTS77' (12) 



so that 



h ~2h s G 

n _Za?P 

2A 3 w' 



(13) 



which gives the desired relation between shear mod- 
ulus and center deflection under a given load. 

3.2. Extension into plastic range 

Equations (1) to (6) apply beyond the elastic range 
as long as plane sections remain plane. Equation 
(8) for the relation between twisting moment and 
shearing stress applies beyond the elastic range also. 
We may replace z in the integral of that equation by 
7, using eq (5) so that eq (8) becomes 



m x 



r ~2~2yl J 



ydy. 



(14) 



The change in load dP corresponding to a small 
change in extreme fiber shearing strain dyi is obtained 
by differentiating this equation in respect to y x with 
the result 



dP 

dy 



2h 2 Cn , , h 2 1 , op , , 2 , 

-=-- -r Tydy+-- 2 T l y 1 =~-(-2P+h 2 T 1 ), 

l 7i Jo 7i 7i 

(15) 

where r x is the extreme fiber shearing stress. Solv- 
ing this equation for t x gives 

The extreme fiber shearing strain can be computed 
simply from measured values of extreme fiber bend- 
ing strain e x or from measured values of deflection 
w a by making use of eq (6) and (12) 



7i=2 €l = 



2hw a 



(17) 
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The stress-strain curve in shear can be derived, there- 
fore, from measurements of Pfe) or P(w a ) by using 
(17) to compute the shearing strain y u plotting 
P(7i), and computing t x from the ordinatcs and 
slope of P(ti) in accordance with eq (16). 

It- is interesting to note that eq (16)[is of Jhe same 
lor 11 1 as the equation 

^ = 6P[^+ 2M ] (18) 

which is derived by Nadai [6] to compute the stress- 
strain curve <ti (eO from the bending moment-strain 
curve M(eO for a rectangular beam of width, b, and 
height, h, in pure bending; t x is the extreme fiber 
bending strain produced by the bending moment, 
M. The identity in form oi eq (16) and (18) might 
be expected because the square plate is bent plasti- 
cally into an anticlastic surface by equal and oppo- 
site bending moments per unit length about axes al 
right angles to each other. 

4. Design of Specimen 

In designing the specimen, consideration must be 
given to practical deviations from the theoretical be- 
havior because the theory neglects edge effects, mem- 
brane stresses, and deformation by transverse shear- 
ing forces. These effects must be considered in some 
detail in order to dimension the specimen so that the 
deviations from ideal behavior can be either neglected 
or taken into account. 

The edge effect will be considered first. In the 
theory, the concentrated forces at the corners are 
replaced by twisting moments uniformly distributed 
along each edge, which are applied by horizontal 
shearing forces rather than the vertical forces shown 
in figure 2. It follows that the theoretical solution 
will hold only for regions of the plate at a sufficient 
distance from the edges, to make Saint-Venant's 
principle applicable. Dimensional reasoning indi- 
cates that the " sufficient distance" should be pro- 
portional to the thickness, h. The factor of propor- 
tionality may be estimated roughly by examining 
cases for which solutions have been obtained from 
the theory of elasticity. 

The plane stress solution for a bar of constant 
width h compressed by a concentrated load at the 
center of the free ends is given in [7]. This shows 
that the maximum compressive stress, which be- 
comes infinite directly under the load, differs less 
than 3 percent from the average stress at a distance 
h from the end. Seewald's plane stress analysis of 
the local stresses near a concentrated transverse load 
acting at the center of a slender beam of rectangular 
section with simple supports at the ends [8] shows 
that these stresses become negligible compared to 
the extreme fiber bending stresses derived from 
simple beam theory at a distance from the load 
greater than the depth of beam. 

This suggests that the errors resulting from the 
application of Saint-Venant's principle will be small 
in regions of the plate at a distance from the edge 
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Figure 4. Distribution of extreme fiber .shear stress 
(Tx'j/Omax along surface of strip with b/h= 10. 

which is greater than h. 'Phis condition is easily 
satisfied for a plate that is thin in relation to its 
lateral dimensions. Unfortunately, it will be shown 
below that a very thin plate is not acceptable because 
it will develop appreciable membrane stresses under 
the loads required to produce yielding in shear. 
The membrane stresses are negligible only if the 
plate has a width not. exceeding a few times the 
thickness. In such a plate, it is probable that the 
deviations from the assumed stress distribution at 
the edges will have an appreciable effect upon the 
magnitude of the stresses in the central portion of 
the specimen. 

A correction for the edge effect in the twisted 
square plate was suggested to t ho authors by their 
associate, Samuel Levy. He pointed out that the 
twisting of an elastic square plate can be considered 
as a special case of torsion of a rectangular plate 
considered by Kelvin and Tait [9]. The elastic 
stress distribution at crosssections of a long rectan- 
gular plate far from the ends is given in [10]. The 
extreme fiber shearing stress is plotted in figure 4 
for the case b/h= 10 which was chosen, after a pre- 
liminary analysis, as a practical compromise for the 
twisted plate. Kelvin and Tait [11] showed that 
the shearing stresses t x * v » parallel to the plane of the 
strip resist only one-half of the external twisting 
moment M t =Pb transmitted by each transverse 
section of the strip. The other half is resisted by 
the transverse shearing stresses r y * t acting on the 
section. 

A redistribution of stresses will take place as the 
strip is twisted beyond the elastic range. The 
nature of this redistribution becomes clear by con- 
sidering the extreme case of stress in a twisted strip 
of perfectly plastic material when all the material 
has yielded. This case was solved by Nadai [12] 
with the sand heap analogy. The resultant stress 
distributions are shown dotted in figures 4 and 5. 
Integration of the stresses for the perfectly plastic 
state shows that, in the plastic as well as the elastic 
range, M t /2 is resisted by the shearing stresses r x * y > 
parallel to the surface of the strip and M t /2 by tin 4 
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Figure 5. Distribution of transverse shear stress (r y ' z ) ma x 
along x'-axis of strip with section b/h — 10. 

shearing stresses r y > z normal to the surface of the 
strip. 

The shear stress distribution of figures 4 and 5 
and the condition that one-half of the twisting 
moment M t is resisted by t x > v > and the other by r rz 
show that the stresses across transverse sections of 
the strip can be approximated b}^ application of 
equal and opposite concentrated forces near the 
ends of the section and a uniformly distributed 
twisting moment along the entire section. 

This distribution of equivalent external forces 
coincides with that shown in figure 2 as acting on the 
edge y' = b/2 of the square twisted plate provided 
that the corner forces P are moved in from the 
corners by a small distance to make P/2 act at the 
same point as the resultant of the transverse shear- 
ing stresses r y > z which are shown in figure 5 for 
points along the x'-axis. It may be concluded that 
the distribution of forces acting on two edges of the 
square plate in figure 2 nearly balances the shearing- 
stresses acting across any transverse section of a 
long strip of width, b, and thickness, h, which is 
twisted by a torque M=Pb applied to the ends. 
The square plate can be regarded, therefore, as a 
portion of such a long twisted strip with a distribu- 
tion of shearing stress across the section y' = which 
approximates that in the transverse sections of the 
strip. Consequently, the extreme fiber shearing- 
stress distributions indicated in figure 4 would also 
apply to the square plate. The shearing stress in 
the central portion would be greater than that com- 
puted from eq (16) because of the above-mentioned 
edge effect. In the elastic case, this increase would 
be about 7 percent for b/h=10, according to figure 4. 

In general, we find from [10] that the extreme 
fiber shearing stress in the central portion of the 
longer sides of a rectangular section can be expressed 

by 



Sk 2 ' 



where k 2 is a factor depending on b/h. 
b/h^>5, the relationship reduces to 



i 



(19) 



For values of 



n 1 -0.63^/6 



(20) 



In the limiting case of a perfectly plastic material 
which has yielded throughout the section, we have 



from the sand-heap analogy [12] 



1 



1-0.333A/6 



For b/h -10, we obtain from (20) and (21) 
n c =1.067ri elastic case} 
t u = 1.0347! plastic case; 



(21) 



(22) 



The actual correction for edge effect beyond the 
elastic range should lie somewhere between the limits 
given by eq (20) and (21). It should be closer to 
(20) as long as the plastic yielding is localized in 
the outer fibers of the plate. 

The effect of membrane stresses will be considered 
next. The stiffening effect of membrane stresses is 
noticeable experimentally during twisting of a thin 
square plate as the transverse displacements w(x,y) 
become comparable to the thickness of the plate. 
Figure 6, taken from Figure 11 of [13] shows the 
effect for a 7.2X7.2X0.126-in. plate of 24S-T 
aluminum alloy. The two straight lines in the 
figure are the strains €i along the diagonals computed 
from eq (9) for the linear elastic case. It is seen that 
the measured load P agrees closely with the calcu- 
lated values at very small strains, but that it deviates 
increasingly with increasing strains. The deviation 
is close to 30 percent for a strain ei =0.0008. 

The effect is analyzed in the appendix for the 
elastic range on the assumption that the transverse 
deflection of the plate continues to be given by 
eq (3). Equation (67) in the appendix shows that 
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Figure 6. Extreme fiber bending strains at center and along 
diagonals of 7.2X7.2X0.126 in. plate of aluminum alloy 
24S-T twisted by transverse forces P applied at the corners. 
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the relative error in extreme fiber bending strain e x 
caused by the stiffening effect, is equal to 






=0.03666 



ay* 



(23) 



The error goes up with the square of the extreme 
fiber strain ei and with the fourth power of the width/ 
thickness ratio, b/h. It may be made small by 
making b/h small. 

Lastly, consideration will be given to the transverse 
shearing stress. The effect of transverse shearing 
stress will be negligible as long as this stress is very 
small compared to the extreme fiber shearing stress 
parallel to the plane of the plate. The average 
transverse shearing stress r r on a section of a cylinder 
at a radius, r, from one of the corners of the plate 
(fig. 2) is 



rh 2 



(24) 



The extreme fiber shearing st icss in the elastic range 
is from (7), (6), and (9) 



Tl = G yi = 2G (l = l/> /r. 



(25) 



The ratio of the average transverse shearing stress 
to the extreme liber shearing stress in the elastic 
range becomes 



t x 37r r r 



(26) 



This requires h/r to be small, that is, the effect 
becomes negligible as the plate is made very thin. 
Equations (20), (23), and (26) were used to esti- 
mate the deviations from ideal behavior caused by 
edge effect, membrane strain, and transverse shear 
for square plates with width/thickness ratios b/h=5, 
10, 20. The results, rounded off to the nearest per- 
cent, are given in table 1. 





Table 1. 








Edge effect (t, h x )—\ 


b/h 


5 


10 


20 


0.13 
.00 
.20 


0.07 
.02 
.10 


0.03 
.38 
.05 


Membrane stiffening effect 

Transverse shear effect r r /r, 


Aei/ei. 







In table 1 (t u /ti) — 1 is taken as the correction 
according to eq (20) for the elastic case. The ex- 
treme fiber shearing strain Yi=2ei is assumed as 0.016 
to place it well beyond the elastic range. The mem- 
brane stiffening effect was computed from eq (23) 
for the elastic range, since the stiffening will be rela- 
tively insensitive to the yielding that takes place 
near the extreme fibers of the twisted plate. The 
average transverse shearing stress was computed for 
a radius r=b/5 from the corner leaving the inner 



portion of the plate 0.2<V/6<0. 7 available for 
strain and displacement measurements. 

Table 1 indicates that a plate with b/h = 10 
provides a fair compromise between the three differ- 
ent disturbing effects that were considered, 

5. Tests 

5.1. Test procedure 

Specimens were made from aluminum alloy 75S-T6 
sheet nominally 0.125 in. thick for which the tensile 
and compressive stress-strain curves had been ob- 
tained both parallel and perpendicular to the direction 
of rolling [14]. This material was selected because of 
its approximately isotropic properties. The speci- 
mens were taken from near the center of the sheet, 
where there was little variation in thickness. The 
specimen selected for trying out the procedure was 
0.1324 in. thick and 1.25 in. square. This gave a- 
width/thickness ratio b/h — dAA, which approximates 
(he value b/h= 10 suggested in the previous section. 
The specimen was large enough to provide space for 
mounting a deflectometer and for attaching strain 
gages at several locations. It had additional mate- 
rial at each corner, as shown in figure 7, to reduce the 
concentration of transverse shearing stress at the 
corners. A conical seat was machined at each corner 
of the square, a pair on one side at the ends of a 
diagonal and a pair on the other side at the ends of 
the other diagonal. 

The specimen S, figure s, was loaded through 
1/8-in. steel balls in the conical seals at the coiners 
with a beam-and-poise screw-type testing machine 
of 50,000-lb capacity used in t ho 5,000-lb range. In 
this manner of loading, the load acted along vertical 
lines through the centers of the balls, which were 
almost a thickness of the sheet from the neutral 
surface. This resulted in the moment becoming 
larger instead of less as the specimen twisted. This 
effect is opposite in sign to the stiffening effect of 
membrane stresses previously mentioned. Two cor- 
ners of the specimen rested upon balls. One of these 
was seated in a block of steel resting on the upper 
head of the testing machine, and the other was seated 
in the block A, figure 8, which was free to move on 




I'/e" 



Figure 7. Location of wire strain gages on plate. 
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legs, D, placed along one diagonal of the plate, while 
their lozenges were connected to a pair of shorter 
legs, E, placed along the other diagonal. One of the 
gages was cemented to the cross piece holding the 
two longer legs. The knife edge of the other gage 
was mounted on a X-in. tube which was inserted in 
the cross piece midway between the two legs. The 
lozenges were actuated by a %-in. rod, which was 
screwed into the cross piece holding the two shorter 
legs and which was reduced at the lower end to pass 
through the tube with ample clearance. The loz- 
enges were set near their lower limits of travel so that 
the vertical component of the clamping force would 
assist in maintaining contact of the two shorter legs 
with the specimen. The deflectometer was mounted 
on the specimen symmetrically with respect to the 
center with the point of one of each pair of legs in a 
small punch mark on each diagonal and the other in 
a scratch along the diagonal. The span was about 
0.7 in. The deflectometer measured the sum of the 
deflections at the center with respect to each pair of 
points or twice such a deflection, since the loading 
was equal and opposite at each pair of corners. With 
this deflectometer, deflections could be read to the 
nearest 0.000002 in. 



Figure 8. Plate under load with deflectometer in position for 
measuring deflections along principal diagonals. 

rollers in the direction of the diagonal. The other 
two corners were loaded by means of a hanger con- 
nected to the loading head of the machine by a rod 
having a ball-and-socket joint at each end. The 
upper cross pieces, B, in which the balls were seated 
are shown in figure 8. They were connected to a 
common piece below the head by steel tapes, C. The 
ball joints at each end of the rod and the pin joints 
at each end of the tapes had the effect of equalizing 
the load on the corners of the plate. 

SR-4 type A-7 strain gages were mounted on the 
specimen in pairs on opposite sides of the specimen. 
The gages were connected with an SR-4 model K 
portable strain indicator through a multiple-pole 
mercury-contact switch. A gage cemented to a 
piece of the sheet material was used as a temperature 
compensating gage. Since each portion of the spec- 
imen was subject to both tensile and compressive 
strains at right angles to each other, the gage-factor 
setting was taken as the axial sensitivity of the gage 
computed from the gage factor provided by the 
manufacturer and the data on axial and transverse 
sensitivities for gages of this type [15]. Strain 
readings were estimated to 0.1 division (0.000001). 
A correction was applied to compensate for the strain 
sensitive wire being farther from the neutral surface 
of the specimen than the surface of the sheet; each 
strain was multiplied by the ratio of the thickness of 
the specimen to the sum of the thicknesses of the 
specimen and the gage. 

A deflectometer using two Tuckerman optical 
strain gages with 0.4-in. lozenges, figure 8, was used 
to measure deflections relative to the center of the 
plate. The Tuckerman strain gages were mounted 
so that their knife edges were connected to a pair of 



5.2. Test Results 

a. Elastic Range 

Wire strain gages were cemented on the specimen 
at locations 1 and 3, and 2 and 4, respectively oppo- 
site them on the other side; figure 7. The speci- 
men and deflectometer were set up as shown in 
figure 8, and the specimen was loaded in the elastic 
range. The value of the shear modulus is from (6), 
(12) 



p Tic Tic T \ c (l 

~ y 1 ~2e 1 ~2w a ti 



(27) 



where ti« is the shearing stress after correction for 
edge effect. Within the elastic range, from (20) and 

(25), 

T ic=T2Ti — ^orrr; ( 28 ) 



or with b/h=9A4, 



'h 2 {\-0.mh!b) 



t u =3.21P/A 2 . 



(29) 



Substitution of t u in (27) gives the shear modulus 
in terms of either measured extreme fiber strain ei, 
or measured center deflection w a . Table 2 shows the 
values thus obtained from the initial linear portion 
of plots of deflectometer readings and wire strain-gage 
readings as a function of load P. 

The last value of G in table 2 was computed from 



G- 



E 



(30) 



2(1+*)' 

by substituting for the Young's modulus an average 
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Table 2. Shear-modulus values measy,red at central portion 
of square plate 



Obtained from— 


G 


Deflectometer. _ 


IblinJ 1 
3.82X106 
3.87 
3.82 
3.84 
3.88 
3.90 


Strain gage 1 

Strain gage 2 

Strain gage 3 

Strain gage 4 

Equation (30) __ 





value E= 10.37 X10 6 lb/in. 2 obtained from tensile and 
compressive tests of the aluminum alloy 75S-T6 
sheet, and for Poisson's ratio v=0.33, the nominal 
value for this alloy. 

Table 2 shows that the experimental values of 
shear modulus agreed within 1.5 percent. The 
average for the wire gages was about 0.7 percent 
greater than the value obtained from the deflection 
data and about 1.2 percent less than the value com- 
puted from eq (30). A lower value from the deflec- 
tion data would be expected because of the effect of 
transverse shear. 

The deflectometer was removed, and gages 5 and 7, 
figure 7, and gages 6 and 8, respectively opposite 
them, were cemented to the specimen. The loading 
was repeated in the elastic range. The results are 
given in table 3. 



Table 3. 



Shear modulus values measured in outer portion 
of square plate 



Obtained 
from gages 


G 


5 
6 

7 
8 


lb/in. 2 
3.53X106 
3.46 
4.04 
3.98 



The results show greater strains at locations 5 and 
6 than at any other location. This indicates a devi- 
ation from the uniform state of stress, which becomes 
serious in the corner portions of the plate. 

b. Plastic range 

The specimen was loaded until a shear strain of 
about 2 percent was reached at the extreme fiber. 
Deflections were measured with the deflectometer 
and strain with the four central wire gages. The 
shear stress was partially corrected for edge effect 
by substitution of n, eq (16), in eq (20). Equation (20) 
was used beyond the elastic range to avoid discon- 
tinuity at the end of the elastic range. It should be 
an adequate approximation in the important region 
near (lie knee of the stress-strain curve because only 
a small proportion of the material would be in the 
plastic range. 

The stress-strain data are plotted in figure 9. 
The values of G determined both from the strain 
data and deflection data were a little over 2 percent 
higher than in the first test, the average value for 
the four wire gages being almost 1 percent above the 
value computed from (30). No explanation could be 
found for this small deviation. 
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Shear strain, r , 



Figure 9. Stress-strain curve in shear, aluminum alloy 
75S-T6 sheet. 

c. Comparison with stress-strain curves under axial load 

It is interesting to compare the stress-strain curve 
in shear with that under axial load. Such a compar- 
ison can be made by deriving the octahedral stress- 
strain curve in each case [16]. The octahedral stress- 
strain curves for both states of stress should coincide 
if the maximum shear strain y on the octahedral 
plane (the plane defined as having a normal which 
makes equal angles with the three axes of principal 
stress) is a single-valued function of the maximum 
shearing stress r on the octahedral plane, regardless 
of the direction and proportion of the principal 
stresses. 

The octahedral stress-strain curve T (y ) was com- 
puted from an average axial stress-strain curve <r(e) 
by substitution in the following formulas (16): 



V2 . r ,-/ a\ 



(31) 



The axial stress a was obtained from tensile and 
compressive stress-strain curves of the material 
parallel to and perpendicular to the direction of 
rolling [14] as the average of the tensile and com- 
pressive stresses corresponding to a given axial 
strain e. The value of G used, 3.93X 10 6 lb/in. 2 , was 
the average determined from the strain data in this 
test. The value of E used, 10.37 X10 6 lb/in. 2 , was 
the average obtained from tensile and compressive 
tests on the sheet. The stress-strain curve in sh (un- 
determined from the strain data, figure 9, was con- 
verted similarly into an octahedral stress-strain 
curve by substitution in 



T Z 



•VI 



TU Yo = 



U 



+Vf( Ti ~S) 



(32) 



The resulting curves of octahedral stress versus 
total strain y and versus permanent strain y ~r /G 
are plotted in figure 10. It will be noted that the 
two stress-strain curves are in good agreement ex- 
cept for deviations up to 5 percent near the knee and 
at shearing strains exceeding 0.009. 

In this connection, mention should be made of 
the much larger discrepancies found in exploratory 
twisting tests of a 0. 125 X 0.625 X0.625-in. plate of 
SAE 1025 steel. These tests were made to extend 
the above comparison to a material with a pro- 
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Octahedral shear strain 



Figure 10. Octahedral stress-strain curve, aluminum alloy 
75S-T6 sheet. 



nounced yield point. The shear modulus was found 
to be in good agreement with the value computed 
from Young's modulus and Poisson's ratio. How- 
ever, the octahedral stress corresponding to the yield 
point obtained from the twisting test was 38 percent 
greater than that obtained from the average of the 
yield points of tensile and compressive specimens of 
the material. Twisting tests of long strips of 
width/thickness ratios equal to 5 and 10 gave 
octahedral stresses corresponding to the yield point 
respectively 17 and 26 percent greater than for the 
average of the tensile and compressive tests. A 
twisting test of a long round specimen 0.59 in. in 
diameter, on the other hand, showed a corresponding 
octahedral stress 15 percent lower. 

6. Conclusion 

The stress-strain curve in shear up to strains of 
the order of 0.01 can be obtained for isotropic flat 
sheet from measurements of deflection or of extreme 
fiber bending strain in the center portion of a twisted 
square plate, provided that a correction is made 
for the edge effect. The shear moduli obtained from 
the twisting test of a 0.13X 1.25X1. 25-in. square 
plate of 75S-T aluminum alloy agreed within 2 
percent with the value calculated from the average 
of the tensile and compressive moduli and an accepted 
value of Poisson's ratio. The octahedral stress de- 
rived from the stress-strain curve in shear for the 
twisted plate agreed within about 3 percent with 
the octahedral stress derived from tensile and com- 
pressive tests of the same material except near the 
knee, where the systematic deviations were as high 
as 5 percent. Much greater deviations were found 
in exploratory tests of mild steel. 

Unfortunately, the method is difficult to perform 
on sheet material because the specimen must be 
small in size. Alternative methods should be inves- 
tigated. An examination of other methods indicates 
that the annulus of constant thickness method can 
be used for determining the stress-strain curve in 
shear, even in the presence of a large radial variation 
in shear stress. Preparations are under way for 
carrying out such tests at the National Bureau of 
Standards. 



7. Appendix 

7.1. Membrane stresses in elastic range 

The membrane stresses in a twisted square plate 
can be estimated from the equations for the state 
of strain in a flat plate with large deflections. These 
are given by Love [17] as 



du , 1 /dwY 
~dx) 



(33) 



du . 1 /c 
_ dv , 1 /dwX 



7xy ~dx + dy + 'dx^y J 



It is assumed that the plate of figure 2 is bent into 
an anticlastic surface so that the transverse deflection 
is given by 

(34) 



w= 



x 2 — y 2 



2R 



just as in the case of small deflections, eq (3). 
The curvature is from (11) 



1 
R z 



26, 

" h ' 



(35) 



The magnitude of the membrane strain depends on 
the displacements u, v parallel to the x, y axes as 
well as on w. 

The displacement u in the ^-direction can be 
approximated by a polynomial of the form 



u=^c i p i y', 

1,3 



(36) 



where i, j are positive integers. The coefficients 
c tj of many of the terms in the series must be equal 
to zero in order to satisfy the following conditions 
of symmetry and antisymmetry in the plate, figure 2: 



(37) 



u(x,y)=u(x,—y) 

u(x,y)=—u(-x,y) 

These conditions are satisfied if i, j are confined to 
i= 1,3,5 . . . 
i=0,2,4 . . . 



(38) 



The displacement v at the point (y, x) must be equal 
to the displacement u at the point (x, y) , so that 






(39) 
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Up to and including terms of t 

-cries (36) and (39) become 



le fifth power, the 

u = c lo x J r c l2 xy 2 + c 3G x 3 + c u nf + c :i2 xhf + c w x 5 



= c w y + c l2 x 2 y + c 30 if + c u xSj + c 32 :r 2 y* + c w y 5 



.(40) 



To satisfy the conditions that the membrane stresses 
<Tx f , Tx'u- be zero along the edge 



b 



(41) 



of the plate, and that the strain energy stored in 
the form of membrane stress must be a minimum, 
the six unknowns Ci , c V2 , c 30 , c Uj c 32 , and c 50 were 
determined as follows. 

Equilibrium of forces on an element of the plate 
bounded by the edge and by planes ?/ — const., 
x= const, shows that the membrane stresses along 
the edge are related to the membrane stresses 
<ti, <r«. r TU by 



_a x -\-a y 



+-. 



_ C X —(Ty i 

rx-v— 2 1 



(42) 



The stresses c x , a y , r xy must be adjusted to make 
c x < and T x > y » in (42) equal to zero. The membrane 
stresses are related to the membrane strains e x , e ?/ , 
y xy by Hooke's Law 

E E \ 

<?*= x ~£_ v2 (€ x +ve y ), a y = 1 _-^ 2 (€ v +ve x )j 



E 



"2(1+,) 



7x V 



f 



(4?) 



where E, pare Young's modulus and Poisson's ratio, 
respectively. 

The membrane strains may be computed from 
the displacements, using eq (33). Carrying out the 
computations and using the following notation for 
convenience 



e=i 



C (1 =C, 



y. 
T 



b 

~~R 



c 2 =c V2 b 2 ; c 3=3c 30 6 2 +^g 2 / 

c 4 =c ll b i ; c 5 =c 50 6 4 ; c,~c 3i b* ) 

lead to 

«*= Co + c 2 r/ 2 + c£ 2 + C iV * + 3c 6 |V + 5C 5 ^ 

<* = Co + c^ 2 + C 3 r, 2 + c 4 £ 4 4- 3c 6 |V + .Serf 
y r!/ =[^ 2 -p 2 +(4c i + 2c,)(e+v 2 )}h ) 



(44) 



(45) 



Substituting (45) in (43) and (43), in (42) and using 
the notation 



1— v 



(4»>) 



gives, for the normal and shearing stresses acting 
on the edges, 

E 



2(1-0 



{2c + (c 2 + C 3 )(e+r, 2 ) + 



(c A + 5c,)(^ + v *) + 6c li e v 2 +a^[4c 2 ~p 2 

(4c 4 +2c e )(e+V 2 )]} 
E 



2(l+v) 



(e-m^-c-2)+(5c 5 - Ci )(e+v 2 )] 



M47) 



The shearing stresses t x - v > will be zero for all points 
£, r\ along the edge of the plate only if the expression 
in the square brackets is equal to zero. This re- 
quires that 



Ci C 2 ) OC5 C4. 



(48) 



The normal stress o x > along the edge can be expressed 
as a polynomial of the fourth order in £ by noting 
from (41) that, for points along the edge of the plate, 



<-vr*- 



(49) 



Collecting terms by powers of i and substituting 
(48) leads to 

--^ ^,=2c +c 3 +|c 5 +-| [-4c 3 -20c 5 + 

a(4r 3 +10c 5 + c 6 -p 2 )] + ^ 2 [4c 3 +30c 5 + 
3c 6 +a(-4r 3 -30c 5 -3c 6 + p 2 )] + 
A 2>[-20c 5 -6c 6 +a(4c 8 + 40c 6 )] + 

r[20c 5 +6c 6 + a(-4c 6 -40c 5 )] = (50) 

The stress v x > will be zero for all points £ along the 
edge only if each one of the five brackets on the 
right is equal to zero. Examination of these brackets 
shows that only three of them are independent of 
each other. One may solve the resulting three 
equations for any three of the coefficients c , c 3 , 
c 5) Cq in terms of the remaining coefficient. 

Solving in terms of c and assuming for Poisson's 
ratio the typical value j> = 1/3 and hence from (46) 
a= 1/2, one obtains from (50) 

2 1 

c 3 =— 3c +p 2 /8; c 5 =-c ——p 2 ; c 6 =0, (51) 

Substitution of (51) and (48) in (45) shows that the 
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strains reduce to 



y xy =[4c 3 -p i + 20c 5 (e + V 2 Mv 



(52) 



It remains to compute the constant c by minimizing 
the membrane strain energy stored in the plate. 
The membrane strain energy per unit volume, 
^0= K< 7 * € *+ <r y e y -{- r xy y X y) is with (43) and (52), assum- 
ing\ = l/3 



V =f^E(8el + yl v ). 



(53) 



The' elastic membrane strain energy V m stored in the 
plate is obtained by integrating (53) over the volume 
of the plate 

(bl-y/Y)- 



Vm = 4 \ n n lohdydx 

Jx = Jy=0 



Substitution of (52) in (53), of the resulting expres- 
sion in (54), and integration leads to 

V m = ^^- [2c§+| c a c 3 +± (c 2 +10c c 5 )+ 

I ^(32c 3 2 -8c3P 2 +P 4 )+ ! ^(36c3 C5 -c 5 p 2 )+i| C l] 

(55) 
The membrane energy V m is minimized by computing 



dc 



-0, 



(56) 



and substituting (51) to obtain an equation having 
Co as the only unknown. Solving for c gives j jj jj 



Co=- I 4gP2=-0.0205p 2 . 



(57) 



The membrane strains and stresses throughout the 
plate are then obtained by substituting (57) in (51) 
and the resulting constants in eq (47) and (52) and 
the latter then in (43). 

The stresses can be plotted in dimensionless form 
by multiplying them by h 2 /b 2 y in, where b/h is the 
width/thickness ratio for the plate and 71, n are the 
extreme fiber shearing strain and stress due to 
bending. These are from (6), (7), and (11) 



Ti r 



h. 



n = G 



(58) 



The results of such computations are plotted in fig- 
ure 11 as solid lines. The plots show that the mem- 
brane stress at the center of the plate consists of a 



uniform compressive stress 



= 0.082 p 71 ti. 



(59) 



The membrane stresses are largest at the center of 
the edge of the plate. The membrane normal 
stresses at this point are 



(r y > = 0.136 -p Tl T l; <T X '=0. 



(60) 



A check on this solution seemed advisable since it 
did not satisfy the conditions of equilibrium at all 
points of the plate; it would be necessary to apply a 
system of body forces per unit volume with compo- 
nents X, Y in the x, y directions in order to maintain 
equilibrium in all plate elements. The equilibrium 
conditions would be satisfied more nearly if c is 
chosen to minimize the average value of X 2 -\-Y 2 for 
the entire plate instead of minimizing the strain 
energy. This leads to 



Co =:_0.0187p 2 , 



(61) 



which is about 10 percent below the previous value. 
The membrane stresses corresponding to this value 
of Co are shown as dotted lines in figure 11. They 
differ less than 10 percent from the stresses previously 
obtained. 

A second check on the adequacy of the theory was 
obtained by computing from it the nonlinear relation 
between load P and extreme fiber bending strain €1 
and comparing it with the experimental values given 
in figure 6. The relation between P and €] for the 
elastic range was obtained by equating the incre- 
ments of work done by the corner loads dW P and the 



h 2 



b 2 yiri 




Figure 11. Elastic membrane stresses in twisted square plate. 
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resulting increase in si rain energy dV consisting of 
the increment </\\, stored in the form of bending 
strains and the increment dV m stored in the form of 
I membrane strains 



J\Y P =dY h + dY 



(62) 



Each one of the three terms in this equation can be 
regarded as a function of the curvature ratio p=b/R. 
It can be shown easily that the work done by the 
corner forces P corresponding to an increment of 
curvature ratio dp is 



dW P =Pbdp. 



(63) 



The increase in bending strain energy can be derived 
from the bending strain, eq (4) as 



dV b 



Eh" 



pdp. 



(64) 



The increment in membrane strain energy is obtained 
from (55) after substituting (57) and (51) 
leads to 

dV m =0.001Ui\ Ehbydp. 



Substituting (63) to (65) in (62), solving 
replacing p by 

p= iT 2 /t 6 " 
gives the desired relation for P{ti), 
„ Eh 2 



This 

(65) 
for /', and 



-[ 



1+0.03666 



m 



(66) 



(67) 



This reduces to eq (9) for small values of e x with 
i/=l/3. Substitution of #=10.4X10 6 lb/in. 2 , h= 
0.126 in., 6=7.2 in. for the plate of figure 6 leads to 
the solid line shown in that figure. The observed 
points come close to this curve, thereby providing a 
second check on the adequacy of this approximate 
theory for computing membrane stresses in the 
elastic range. 
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